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Abstract 

The ejection of DNA from a bacterial virus ("phage") into its host cell is a bio- 
logically important example of the translocation of a macromolecular chain along its 
length through a membrane. The simplest mechanism for this motion is diffusion, 
but in the case of phage ejection a significant driving force derives from the high de- 
gree of stress to which the DNA is subjected in the viral capsid. The translocation 
is further sped up by the ratcheting and entropic forces associated with proteins 
that bind to the viral DNA in the host cell cytoplasm. We formulate a general- 
ized diffusion equation that includes these various pushing and pulling effects and 
make estimates of the corresponding speed-ups in the overall translocation process. 
Stress in the capsid is the dominant factor throughout early ejection, with the pull 
due to binding particles taking over at later stages. Confinement effects are also 
investigated, in the case where the phage injects its DNA into a volume comparable 
to the capsid size. Our results suggest a series of in vitro experiments involving the 
ejection of DNA into vesicles filled with varying amounts of binding proteins from 
phage whose state of stress is controlled by ambient salt conditions or by tuning 
genome length. 



1 Introduction 

A crucial first step in the life cycle of most bacterial viruses involves binding of the virion to 
a receptor protein in the host cell membrane followed by injection of viral DNA. The viral 
genome is typically about 10 microns long, and its translocation from outside to inside the 
host cell is accomplished over times that vary from seconds to minutes. The wide range of 
mechanisms responsible for injection of phage genomes has recently been systematically 
reviewed including many references to the last few decades of relevant literature. 

In the present paper we formulate a general theory of chain translocation that takes 



into account many of the physical processes involved in actual phage life cycles. These 
processes include: diffusion of the DNA chain along its length; driving forces due to stress 
on the DNA inside the viral capsid; resisting forces associated with osmotic pressure in 
the host cell; cell confinement effects that constrain the injected chain; and ratcheting 
and pulling forces associated with DNA-binding proteins in the host cell cytoplasm. 

Considerable effort has been focused on the energetics of packaging and ejecting DNA 
in phage. In particular, theoretical work 0, H, S S B S HO, III has shown that the dom- 
inant source of stress on the DNA in the capsid results from strong repulsive interactions 
between neighboring portions of double helix that are confined at average interaxial spac- 
ings as small as 2.5 nm. Another major contribution comes from the bending stress that 
arises from the capsid radius being smaller than the DNA persistence length. The force 
needed to package the genome against this resistance is provided by a virally encoded 
motor protein that pushes in the DNA along its length. Recent laser tweezer measure- 
ments have confirmed that this force increases progressively as packaging proceeds, 
i.e., as the chain becomes more crowded and bent, reaching values as large as 50 pN 
upon completion. Conversely, the force ejecting the DNA upon binding of the phage to 
its membrane receptor has been shown [Ij, to decrease monotonically from tens of 
piconewtons to zero as crowding and bending stress are progressively relieved. In the 
current paper we consider the dynamics of phage ejection and attempt to distinguish the 
relative importance of these large, varying, "internal" forces and the binding particles in 
the external solution (bacterial cytoplasm). 

It is useful at the outset to consider the simple diffusion limit of the translocation 
process. More explicitly, consider the case in which a chain is threaded through a hole in 
a membrane dividing one solution from another. If the chain is free, i.e., in the absence of 
pushing or pulling forces and of binding particles, it will simply diffuse along its length, 
experiencing a friction associated with its passage through the membrane and the viscosity 
of the solution. The time required for its translocation from, say, the left to the right 
will be L'^/2D = r^, where L is the length of the chain and D is its effective translational 
diffusion coefficient. 

Suppose now that particles are added to the right-hand solution which bind irreversibly 
to the chain at regularly spaced sites as soon as they diffuse into the solution. Then, if s 
is the spacing between these binding sites, the diffusion of the chain will be ratcheted each 
time another length s has entered the solution corresponding to the fact that 

the chain cannot move backwards through the hole at a site where a particle is bound. 
Accordingly, the time it takes for the entire chain to appear on the right is simply given by 
s'^/2D - the time required for diffusion between a pair of neighboring binding sites - times 
the total number of sites, L/s. It follows that the overall translocation time in the presence 
of perfect ratcheting is reduced by a factor oi s/L over that for free diffusion. When the 
binding of particles is reversible - they do not remain bound indefinitely, thereby allowing 
some sites to diffuse backwards through the hole - the translocation time is increased by 
a factor of (1 + 2K) compared to p erfect ratcheting, where K is the ratio of "off" and 
"on" rates for particle binding |l5l Il6i|. Finally, note that the ideal ratcheting time of 
Ls/2D corresponds to a velocity of 2D / s and hence, by the Stokes-Einstein relation, to 
a force of 2kBT / s pulling the chain into the particle-containing solution [l7i |. 

When the particle binding is reversible, however, it turns out that there can be an 



additional correction to the ratcheting dynamics, one that can significantly shorten the 
translocation time below Ls/2D = rideairatchet- This effect requires that the dffusive motion 
of the chain is slow enough and is due to the fact that the entropy of reversibly bound 
particles increases when there is more chain for them to explore. As a result, the entropy is 
an increasing function of chain length available in the right-hand solution. Indeed, in the 
limit of fully equilibrated binding, the system is equivalent to a one-dimensional Langmuir 



adsorption problem [17|,|18|(P. G. de Gennes, personal communication) . More explicitly, 
the ID Langmuir pressure can be written in the form Pio = (fc^T/s) ln{l + exp[(e -|- 
fi)/kBT]}, where e > is the energy lowering of the adsorbing particles upon binding and 
fi is their chemical potential in solution. Note that in the limit of large binding energy 
((e + jj)/kBT ^ 1) this pressure reduces simply to (e + fi)/s, which - because pressure is 
force in a ID system - can be directly interpreted as the force pulling on the chain due 
to the reversible binding of particles. Note further, in the large binding energy limit, that 



this force is necessarily large c omp ared to the ideal ratcheting force, 2A;^T/s [Jj 

Ambjornsson and Metzler (l8| have recently clarified the various timescales that de- 
termine the different regimes of chain translocation in the presence of "chaperones" , i.e., 
binding particles. The first, tq, is the time needed for the chain to diffuse a distance of 
order s, the separation between binding sites. The second and third are Tocc and Tunocc, the 
characteristic times that a binding site remains occupied and unoccupied, respectively. 
Tocc and Tunocc are related by the equilibrium relation. 



exp 



e + fi 



''"unocc \ 

Finally, Tunocc cau be approximated by the typical time it takes for a particle to diffuse a 
distance of order i? (~ Cq ) between binding free particles: 

r = ^ — (2) 

where Dq is the diffusion coefficient of the particles. One can then distinguish between 
three different regimes: 

1. Diffusive regime: tq <C Tunocc, Tocc- Here the binding particles are irrelevant to the 
chain translocation because the chain diffuses its full length in a time too short for 
the particles to bind. 

2. Irreversible binding regime: Tunocc <^ Tq <^ Tqcc- Here particles bind essentially 
irreversibly on a time scale short compared to the time it takes for the chain to 
diffuse a distance between binding sites. We shall refer to this as the "ratcheting" 
regime. 

3. Reversible binding regime: Tunocc, Tqcc ^ tq. Here diffusion of the chain along its 
length is slow compared to the time required for an "on" /"off" equilibrium of the 
binding particles to be achieved. We shall refer to this as the "Langmuir" regime. 



It is also important to clarify some relevant length scales involved in the problem. 



Specifically, we distinguish between two extremes of how the separation, s, between bind- 
ing sites compares with the range, 6, of the attractive interaction between binding particle 
and chain. Pure and "perfect" ratcheting will arise when Tunocc <^ tq <C Tqcc, independent 
of the relative values of S and s. "Imperfect" ratcheting will arise when Tunocc, Tocc <^ To, 
but 6 <^ s. The translocation time for the "imperfect" ratchet is higher than the "perfect" 
ratchet by a factor of (1 + 2K). Finally, when Tunocc, Tqcc ^ tq and 5 ^ s, in addition 
to the "imperfect" ratchet we also have a Langmuir force on the chain. Note that if the 
binding free energy between DNA and the binding proteins is very large then K <^1, and 
the imperfect ratchet is no different than the perfect one. In this paper we will always 
take this limit. 

Before proceeding further it is instructive to make some numerical estimates. Within 
this simple translocation model all time scales are naturally referenced to that for pure 
translational diffusion of a chain along its length, and hence to the diffusion coefficient 
D introduced earlier. In reality, however, the DNA ejection process is enormously more 
complicated, since the chain moving through the tail of the phage is feeling not only 
the friction associated with the few hydration layers surrounding it but also the viscous 
effects arising from interaction with the inner surface of the tail just nanometers away. 
Furthermore, this chain portion is connected to the lengths of chain inside the capsid and 
outside in the cell cytoplasm. The chain remaining inside the capsid moves by reptating 
through neighboring portions of still-packaged chain and/or by overall rotation of the 
packaged chain. All of these latter motions involve viscous dissipation that is insufficiently 
well-characterized to enable realistic estimates of diffusion time scales, even though one 
can distin guis h between different dependence on chain length for each of these dynamical 
processes |l9|,|9|. As a result, in pursuing the simple translocation picture as a model for 
overall phage ejection kinetics, we resort to using an effective diffusion coefficient D to 
define the unit of time, = L^/2D. 

A strong upper bound for D can be obtained by considering the part of the dissipation 
arising as the chain moves through the tail portion of the virus. Taking into account only 
the friction between the DNA and the fluid in the tail we have, for example [i^ Il9l |. 
( = 27rlr]/ \n{A/d). Here ( is the friction coefficient, / is the length of the tail, rj is 
the viscosity of water, A is the inner diameter of the tail, and d is the diameter of 
the double-stranded DNA. Taking / = lOOnm, rj = 10"^pN-s/nm^, A = 4nm and 
d = 2nm, we find ( = 9 x 10~''pN-s/nm and hence a diffusion coefficient {D = kBT/Q 
of 5 X lO^nm^/s. For a typical phage genome length (L) of 10/im, this in turn leads to 
a diffusional translocation time (r^ = L'^ /2D) of about 10 seconds, not unlike ejection 
times measured for phage A 0. Recall, however, that this estimate .s based on a value 
for D which is a strong upper bound, because of all the viscous dissipation contributions 
that were neglected, suggesting that the actual unassisted diffusional time is likely several 
orders of magnitude larger than this 10 seconds estimate. Indeed, the outcome of the 
work presented below is that the translocation time is shortened beyond by several 
orders of magnitude by a combination of effects dominated by pressure in the capsid and 
binding particles in the external solution. A schematic of the role of these various effects 
is shown in Fig. ^ 

The outline of our paper is as follows. In the next section we include the effect of 
capsid pressure by formulating a Fokker-Planck description of translocation driven by a 
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Figure 1: Schematic showing the various physical effects which assist bare diffusion in 
the process of phage DNA ejection. The DNA cross-section is not shown to scale: its 
diameter is 2 — 4 nm, as compared with a capsid interior diameter that is ten times 
larger. The spring denotes schematically the stored energy density resulting in a force F 
acting along the length L — x oi chain remaining in the capsid. The small spheres denote 
particles giving rise to an "external" (cytoplasmic) osmotic pressure Ilosmotic, while the 
green particles labeled i and i + 1 are successive binding particles. 



combination of diffusion and spatially varying force, i.e., a force pushing the chain from 
one side that depends on the length of chain remaining on that side (corresponding to the 
portion still in the capsid and hence experiencing stress due to crowding and bending). 
We evaluate the mean-first-passage-time (MFPT) for translocation of an arbitrary length 
and thereby calculate the length ejected as a function of time, using estimates of the 
spatially-varying ejection force from recent theories of phage packaging energetics. We 
find that the translocation times are 2 to 3 orders of magnitude faster than the diffusional 
time. We also treat the case of ejection into a volume comparable to the capsid size 
(mimicking, say, studies in which phage are made to eject into small vesicles that have 
been reconstituted with receptor protein j23l.l23|) and find that the ejection time decreases 
from its maximum value when the confinement scale is made larger or smaller than the 
capsid size. In SecElwe treat the further speed-up in translocation due to ideal ratcheting 
and to reversible particle binding, respectively. We find that the simple ratcheting effect 
is small compared to that arising from the entropic force of reversible particle binding. 
The effect of reversible particle binding decreases the translocation time by another order 
of magnitude beyond that due to capsid pressure effects. We conclude in Sec. El with a 
discussion of related work by others, of additional contributions to ejection dynamics that 
will be studied in future theoretical work (in particular, the effect of RNA polymerase 
acting on the ejected DNA), and of experiments planned to test the various predictions 
made in the present work. 



2 Kinetics of ejection driven by packaging force 

As discussed in Sec. we focus here on a chain which has been confined in a viral capsid 
and which is ejected from it through a hollow tail just big enough to accommodate its 
diameter. To elucidate the essentials of this ejection process we describe the translocation 



of the chain as a "diffusion- in-a- field" problem 2^ 2^ 2^. In the present case, involving 
the translocation of a linear polymer along its length, the diffusion coordinate is a scalar, 
i.e., the length of chain x that has been ejected from the tail of the virus. The external field 
is described by the potential energy U{x) that gives rise to the force F{x) = —dU{x)/dx, 
pushing on the chain when a length x of it has been ejected. This force is due to the 
remaining chain length L — x being confined inside the capsid and thereby subjected 
to strong self-repulsion (f/rcp) and bending (?7bend)- The corresponding potential U{x) = 
Urcp{L — x) + Uhwd[L — x) is the free energy calculated in recent theories of DNA packaging 
in viral capsids [8|, |7|, . This energy is seen to decrease dramatically as ejection proceeds 
(i.e., as X increases), and so does the magnitude of its slope that constitutes the driving 
force for ejection. 

The one-dimensional dynamics of a diffusing particle in the presence of an external 
field is a classic problem in stochastic processes j28[, and, as argued above, can be tailored 
to treat the translocation of phage DNA under the action of an ejection force F{x) = 
—dU{x)/dx. Accordingly, the probability p{x,t) of finding a length x ejected at time t is 
given by the Fokker-Planck equation 

dp{x^_ d ( dpix,t) D dUix) \ 

As part of this stochastic description of the translocation-under-a-force process, it is 
natural to define a mean-first-passage-time (MFPT), that gives the average time it 
takes for a length x to be ejected in the presence of the external field U{x), namely jl^. 



t{x) = —J dxi exp 1^ — J J exp 



(4) 



It is useful to consider several limits of this general equation, the first corresponding to 
the familiar case of no external field. From ?7 = the integrals in MFPT reduce trivially 
to x^/2D, giving the expected diffusion time, t{x) = x'^/2D. 

For the case of constant force, i.e.^_^ = —Fx + constant, the integrals in MFPT can 
also be evaluated analytically, giving [l6[ 

, . _ x^ exp[-pFx] +pFx-l 

^ConstantForcc (.-^j ^ (/Ji^a;)^ ' 

Here we have written /3 for l/ksT, and taken F = —dU{x)/dx > to denote the constant 
force driving translocation of the chain to the right. In Sec. El we will apply Eq. El locally, 
over each segment of length s associated with a binding site, to calculate the ideal ratchet- 
ing corrections to force-driven translocation. Note that simple and ratcheted diffusion are 
overwhelmed by force-driven translocation when f3FL ^ 1 and f3Fs ^ 1, respectively. 

In the most general instance of spatially varying "external" field U{x), as in the 
case of capsid-pressure-driven translocation, the integrals in Eq. |31 must be evaluated 
numerically. In this way we calculate t{x) from Eq. |31 for the U (x) determined from a 
recent treatment 0, EH of the packaging energetics in phage capsids. This provides a 
one-to-one correspondence between each successive time t{x) and the fraction of chain 



ejected x(t)/L at that instant. 

In Purohit et al. 0, 0| the DNA inside the phage capsid is assumed to be organized in 
a hexagonally packed inverse-spool. The potential U{x) is expressed as a combination of 
the bending energy and the repulsive interaction between the DNA strands, and is given 
by, 

U{x) = Urcp{L - X) + Ul,cnd{L - X) 

= V3Fq{L-x){c^ + cd)exp{-d/c) 
Vid Jr,„ r 

Fq and c are experimentally determined constants describing the interaction between 
neighboring DNA strands, ^ is the persistence length of DNA, d is the inter-strand spac- 
ing, i?out and -Rin are the radius of the capsid and the inner radius of the DNA spool, 
respectively, and N{r) is the number of hoops of DNA at a distance r from the spool 
axis. We are interested in finding the internal force on the phage genome as a function of 
genome length inside the capsid. We do so using Eq. Eland simple geometrical constraints 
on the phage genome inside the capsid. The number of loops N{r) in Eq. IHl is given by 
z{r)/d, where z{r) = {R'^ut " r^')^!'^ is the height of the capsid at distance r from the 
central axis of the DNA spool. The actual volume available for the DNA - V^(i?in, -Rout) ^ 
can be related to the genome length L — x in the capsid, and the inter-strand spacing c?, 
giving an expression for R\a in terms of c?, i?out and L — x. This relation can be substituted 
for i?in in Eq. IHl which then can be minimized with respect to d to give the equilibrium 
inter-strand spacing as a function of the genome length L — x inside the capsid. In this 
way we determine the total packing energy as a function of genome length inside the 
capsid (L — x) or as a function of the DNA length ejected x, i.e., U{x). Using this result 
and Eq. |3]we can evaluate the MFPT, for the DNA ejection in A as a function of 
the length ejected. The corresponding fraction ejected, x(t)/L, is shown as a function of 
time in Fig. |21 with the label "no confinement" . 

The value of D can be estimated on the basis of this simple model by the following 
procedure. The experiment by Novick and Baldeschwieler 2^1 showed that in a buffer 



containing 10 mM of Mg"*"^ it took roughly 50 seconds for phage A to completely eject its 
genome. The values for Fq and c in buffers containing Mg^"*" have been measured j3ol |. 
Since the values measured for 5mM and 25mM Mg^"*" were not significantly different, we 
assume that the forces at lOmM will be identical, i.e., Fq = 12000 pN/nm^ and c = 0.3 nm. 
Using these values in Eq. EJand numerically evaluating it for x = L = 48500 x 0.34 nm 
we find the total time for A to eject its genome of 48.5 kbp is t ~ (lO^nmVD)seconds. 
Then, since this value is experimentally estimated to be around 50 seconds [2j|, we infer 
that D ^ lO'^nm^/s. This is about 3 orders of magnitude smaller than the D estimated 
in Sec. consistent with all the sources of dissipation that were left out of that estimate. 

An interesting application of our estimates is to the experiments in which viruses eject 
their DNA into lipid vesicles 0, 24, |^. Here lipid vesicles are reconstituted with 



the receptors recognized by the phage of interest, and then mixed with a solution of the 
phage. The phage binds to the receptor and ejects its DNA into the vesicle. We argue 
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Figure 2: Ejection time for phage-A injecting its genome into vesicles of radius 29, 50 and 
100 nm. The capsid radius of the phage is 29 nm. It can be seen that the amount of DNA 
injection increases as the ratio of the vesicle radius to the capsid radius increases. On the 
time scale depicted here, there will be essentially no ejection due to pure diffusion (which 
takes place instead at times of order 1, in units of L^/D). 



that the amount of DNA ejected into the vesicle and the corresponding time depend on 
the radius of the vesicle. In particular if the vesicle has a radius comparable to that of 
the viral capsid there will be a build-up of pressure inside the vesicle due to the ejected 
DNA. Ultimately, the ejection process will come to a halt when the force on the DNA 
from the capsid equals the force from the vesicle side - this can be thought of similarly 
from the free energy perspective as a free-energy minimizing configuration. Hence, the 
ejection will not, in general, be complete. 

We can work out the ejection rate for this process as follows. If x is the length of 
genome ejected into the vesicle, we denote the free energies of the DNA inside the viral 
capsid and the vesicle by UcapsidiL — x) and U^esicieix), respectively. The total free energy 
will be given by, 

U{x) = f/capsid(-^^ -X) + f/vcsicle(x). (7) 

As explained before, we already know f/capsid(-^ — x) - see Eq. IHl the expression for 
fA'csicio (a^) can be obtained similarly by assuming that the vesicle is like a spherical capsid 
and the DNA configuration inside is similar to that inside the viral capsid. Our assumed 
structure for the DNA in the vesicle is a highly idealized model, though we note that 
electron microscopy on such vesicles demonstrates that DNA within them can adapt 
highly ordered configurations ji^. In the limit where the vesicle radius is large compared 
to that of the phage capsid we will recover the free injection result (DNA ejecting from 
phage into the surrounding solution). 

The injection process will stop when the total free energy reaches a minimum, i.e., the 



total force on the DNA is zero. The time for DNA injection is given by Eq. HJ We have 
worked out the kinetics of the ejection for the bacteriophage A (radius ~ 29 nm) ejecting 
its genome into vesicles of radius 29, 50, 100 nm. The phage is taken to be suspended in 
a solution of Mg"*"^ ions, and similarly the vesicle, with concentration that approximately 
gives the same values for Fq and c as discussed earlier. This yields a prediction for the 
kinetics of injection for different vesicle radius. It can be seen from Fig. |21 that when 
the size of the vesicle is comparable to the capsid size there is only a partial ejection of 
the DNA. When the vesicle size is almost twice the size of the capsid nearly the entire 
genome is ejected, except for the last part of the DNA, which takes "extra" time because 
of the resistance offered to it from the DNA inside the vesicle. Finally, when the vesicle 
is more than three times the size of the capsid, DNA gets completely ejected from the 
phage capsid as if there were no vesicle. It is interesting to note that in the initial stages 
of ejection all the curves for various vesicle sizes fall on one another because there is no 
resistance to the injection, but as the ejection proceeds each curve reflects a different 
resistance. 

It is also possible that the arguments given above for in vitro ejection into vesicles 
could be relevant to thinking about ejection into the crowded environment of a bacterial 
cell I33I 1 . As a result of the crowding within the host bacterium, the viral DNA may 
be subject to confinement effects like those induced by vesicles. 

3 DNA ejection in the presence of DNA binding pro- 
teins 

The E. coli cell has as many as 250 types of DNA binding proteins js^l • Some fraction of 
these proteins likely bind either specifically or non-specifically to the phage genome as it 
enters the host bacterium. Accordingly, we consider what happens if the phage DNA is 
swarmed with binding proteins upon its entry into the host cell. Depending on the binding 
on/off rates, binding site density, and the strength of binding, we have a corresponding 
speed-up of the DNA injection into the bacterial cell. In this section we explore this effect 
and see how, in addition to the speed-up, it helps the phage inject its DNA against the 
osmotic pressure in the host cell. 

Throughout the following analysis of particle binding effects, we assume that the chain 
is stiff on length scales (e.g. lO's of nanometers for double-stranded DNA genomes) large 
compared to the size of the relevant binding particles (typically a few nanometers). We 
also assume that the binding particles are comparable in size to the distance between sites; 
for an estimate of Langmuir forces in the more general case of larger binding particles, 
see 0. 

3.1 DNA ejection due to the ratchet action. 

Consider a scenario in which host cell binding proteins irreversibly bind on to the DNA at 
a rate much faster than the translocation rate. In such a case, once a binding site is inside 
the cell, it is immediately occupied by a binding protein. If the protein stays bound long 
enough, compared to the translocation time, it will prevent thermal fluctuations from 



retracting the DNA back into the capsid. As a result, the DNA will diffuse only between 
consecutive binding sites, instead of along its complete length. Depending on the spacing 
between the consecutive sites, it will bring about a speed-up in the translocation compared 
to when it is only force-driven [l6| . 

For simplicity we assume that the protein binding sites are uniformly distributed along 
the length of the genome. If the distance between the consecutive binding sites, s, is small 
compared to genome length, i.e., L ^ s, we can assume that the internal force on the 
genome due to the packaged DNA is effectively constant while the DNA chain is diffusing 
between sites. In that case the MFPT, ti, for the DNA to translocate the distance s 
between the binding sites i — 1 and i is simply given by Eq. with x replaced by s, and 
F replaced by Fj. The internal force F is of course a varying function of ejected length 
X, but to a good approximation is constant over each interval of length s. The subscript 
i on the force F denotes this approximately constant force on the DNA chain when the 
translocation is taking place between the i — 1 and i^^ binding sites, i.e., when length 
(z — l)s has been ejected. 

The total translocation time for ejecting length x of the DNA is given by a sum over 
the MFPTs for all the sections of length s, along the length x ejected. The MFPT as a 
function of x is given by. 



The corresponding plot for the fraction ejected, x{t)/L, as a function of time is shown in 
Fig. 121 for s = 20 nm: the ratcheting reduces the injection time by half as compared to 
when the ejection results exclusively from the internal force. From Eq. it can be seen 
that the time will decrease exponentially as the spacing s decreases. But, since s can not 
be smaller than the size of the binding proteins, which is of the order of 10 nm, we have 
a lower bound on s. The important qualitative consequence of the ratchet will be seen 
(see Sec. 13. 3p to be its helping with internalization of the complete phage genome against 
osmotic pressure, when internal force alone is insufficient to carry it out. 

3.2 Reversible force from the binding proteins 

Consider another extreme scenario where DNA injects into a reservoir of binding particles 
and the rate of translocation is slow compared to the time required for the particles to 
bind and unbind from the DNA. In this case, the binding proteins will come to equilibrium 
with the DNA. As a result there will be an adsorption force pulling on the DNA, given 
by |i3l(P- G. de Gennes, personal communication). 
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Figure 3: The fraction of DNA injected in phage A as a function of time (in units of 
L"^ / D) in the presence of binding particles that form a ratchet. The DNA injection purely 
due to the internal force is used as a benchmark, and the spacing between the binding 
sites s = 20 nm. It can be seen that the ratchet reduces the translocation time. The time 
required to internalize the genome solely by the ratcheting mechanism (see lower, straight 
line) is around twice the time taken for the purely internal force-driven mechanism. 



Here /i is the chemical potential maintained by the reservoir of binding proteins, e is the 
binding energy of the proteins with the DNA, and e + yU = AG'(> 0) is the binding free 
energy for the proteins. This adsorption force is the ID Langmuir pressure discussed in 
Sec. H Now, since we have assumed that AG ^ ksT , the binding proteins are mostly 
bound "on" the DNA and very rarely "off" . So in addition to the force there will also be 
a Brownian ratchet as discussed earlier. In order to evaluate the MFPT we follow exactly 
the same process as in the previous section with the addition of this Langmuir force AG/ s 
to Fi. The total MFPT is then given by Eq. Hwith Fi replaced by Fi + AG/s. We take 
a typical value of non-specific DNA-protein binding free-energy of AG = ^ksT js^. The 
plot corresponding to s = 20 nm is shown in Fig. |3 It can be seen that the Langmuir 
force speeds up the genome translocation by almost an order of magnitude. Not only 
that, but even if we do not have an internal force, this mechanism (see Pure Langmuir) 
will inject the complete genome faster than the internal force-driven mechanism. This is 
because after about 50% ejection, the internal force begins to drop below the constant 
value of the Langmuir force. Indeed, from Fig. |31 we see that it is at an ejected fraction 
of about 0.5 that the slope of the "Internal Force" curve drops below the constant slope 
(rate) of the "Pure Langmuir" plot. 

The two cases we described are really two opposite extreme cases for the DNA binding 
proteins. In reality the rate of binding and the equilibration times may not be very fast 
(compared to translocation times) and the translocation rates would lie somewhere in 
between the rates evaluated in this section; for these cases it is necessary to treat the 
dynamical coupling between particle binding and chain diffusion [l7i |. 




Fig. 4 Time in units of L^/D 

Figure 4: The fraction of DNA injected in phage A in the presence of binding proteins that 
bind reversibly, as a function of time (in units of L^/ D.) The presence of reversible binding 
proteins result in a pulling Langmuir force (see text). This pulling force significantly 
enhances the DNA ejection rate over that of the purely force-driven mechanism, by almost 
a factor of ten. 



3.3 Binding proteins enable DNA ejection against osmotic pres- 
sure 

Due to macromolecular crowding jl^l, the E. coli bacterium has internal osmotic pres- 
sures of about 3 atm j^. The work of Evilevitch et al. showed for phage A that 
the ejection process can be partially/completely inhibited by an application of osmotic 
pressure. Hence, it appears that if the phage were to rely entirely on the driving force 
due to the packaged DNA to eject its genome, the time scale for full ejection would be 
prohibitively long. On the other hand, since we know that the genome is completely in- 
ternalized it seems likely that the particle-binding mechanisms described above may play 
a key role in in vivo DNA translocation. In this section we will see that the task can be 
accomplished by the Brownian ratchet and the ID Langmuir force mechanism discussed 
in the preceding Sees. 13.11 & 13.21 

To see how the Brownian ratchet can internalize the genome against the osmotic 
pressure, we use the following procedure. If the osmotic pressure in the host cell is 
Hosmotic, the resisting force acting on the DNA can be approximated |ll|, |8| by 

-^osmotic 

nosmoticTT-RoNA' wherc -Rdna IS the radius of the DNA (about Inm). For an osmotic 
pressure of 3 atm the osmotic force is then estimated to be around 1 pN. We can now 
replace the term F in Eq. El with Fj — -Fosmotic to evaluate the MFPT, tj for the injection 
of the DNA segment between binding sites i — 1 and i. This time ti is then summed over 
all z, as in Eq. |Hl to give the time t{x) and hence x{t)/L. This fraction is plotted in Fig. El 
for the case of spacing s = 20 nm between binding sites, and for an osmotic pressure of 3 
atm. 

It can be seen from the figure (bottom curve) that the time required for internalizing 
the genome is comparable to the time it takes for phage to inject its genome purely by the 



internal force, when there is no osmotic pressure. The internal force for A (data not shown) 
at around 50% DNA ejection is approximately 1 pN, i.e., of the order of Fosmotic- It can be 
seen from Fig. El that the slope of the curve showing ejection in the presence of ratcheting 
and osmotic pressure starts decreasing at that percentage of ejection. The average force 



produced by a Brownian ratchet is 2kBT/s ^ 0.4pN for s = 20nm [17|, |16|. At 60 — 70% 
ejection the internal force is around 0.5pN; the total driving force is then approximately 
0.5 + 0.4 = 0.9 pN, which is almost the same as Fosmotic- This force hence works to eject 
the genome against the external osmotic force. When around 15% of the genome is left 
in the phage capsid, the internal force is almost zero. At this point there is only a small 
amount of the genome still to be ejected and a small differential of -Fosmotic — -^ratchet ~ 0.5 
pN to be worked against. This is accomplished by the Brownian motion of the DNA. 

Now take the second case, when in addition to forming a ratchet we have a ID 
Langmuir pressure, as discussed in Sec. 13.21 To include the effect of the osmotic pres- 
sure we have to subtract the osmotic force -Fosmotic = nosmoticTT-RDNA from the driving 
force Fi + AG/s and use the result in Eq. El This means that, so long as {AG/s — 
nosmoticTT-RoNA) — 0' ^^^^ always have DNA ejection faster than or the same as that 
for the purely force driven non-osmotic pressure case. For the numbers we took in the 
preceding sections, AG/s = 8kbT/s ~ l.GpN, which is greater than -Fosmotic ~ IpN. This 
implies that the phage would inject its genome faster than in the purely pressure-driven 
mechanism. 




Fig. 5 Time in units of LVD 



Figure 5: The fraction of DNA injected in phage A as a function of time (in units of L^/ D) 
for the case in which there is a resistive force due to osmotic pressure. We compare the 
roles of the Langmuir force and the ratchet effect in ejecting the phage DNA against 
osmotic pressure. The spacing s is taken to be 20nm and the osmotic pressure in the cell 
is around 3atm. It can be seen that the Langmuir force easily pulls the DNA against this 
pressure. The DNA translocation by the Brownian ratchet requires a much longer time, 
but it still suceeds in pulling out the genome at time-scales not much longer than the 
ejection by internal force alone with zero osmotic pressure. 



4 Discussion and Conclusion 



This paper addresses the problem of the kinetics of phage injection and the various mech- 
anisms responsible for it. We make use of the available experimental data, existing mod- 
els for phage packaging, and classical Fokker-Planck theory, to make predictions about 
translocation rates for phage DNA ejection. The key quantitative predictions described 
in this paper are: 

• Dependence of ejection rates on driving pressure. As shown in Fig. |2l the driving 
force due to the packaged DNA speeds up the ejection process by 2 — 3 orders of 
magnitude over free diffusion, and thus is a major contributor to the process of 
injection. Also, in the in vitro setting, the smaller the vesicle into which ejection 
occurs, the smaller the amount of DNA injected. In addition, for genomes of the 
same size, the time required for the ejection of the DNA is larger than when into a 
bigger vesicle. 

• Dependence of ejection rates on the presence of irreversible DNA-binding proteins. 
Ratcheting enhances the DNA ejection rate from the viral capsid. The speed-up 
is minor when compared to internal force-driven ejection (see Fig. E)), but as seen 
from Fig. El it is sufficient to pull out the genome against osmotic pressures of up to 
3 atm found inside the bacterial cell. 

• Dependence of ejection rates on the presence of reversible binding proteins. The 
reversible binding of proteins exerts a ID Langmuir force on the DNA. It can be 
seen from Fig. |3] that the presence of this phenomenon significantly enhances the 
DNA ejection rate beyond that due to pressure in the viral capsid. From Fig. El it is 
clear that this force is sufficient to efficiently internalize the phage genome against 
osmotic pressures of up to 3 atm in the bacterium. 

We have several biological examples in mind when we treat these ejection mechanisms. 
In bacteriophage T5 the DNA injection occurs in two steps. The first step transfer, 
which involves ejection of around 10% of the phage genome, is driven by the internal 
force There is then a brief pause, when a protein is synthesized that is implicated 
in the degradation of the host chromosome, thereby freeing the large number of proteins 
that had been bound to it. These latter proteins are now available for binding to the 
injected portion of the phage genome and for pulling the remaining DNA into the cell, 
via the ratcheting and Langmuir mechanisms. 

Similar ideas to those proposed here might also prove useful in those cases where the 
viral genome is translocated as a result of the binding of motor proteins which themselves 
translocate along the DNA. One such example is the pulling force by the NTP-driven 
RNA polymerase (RNAP). RNAP is a very strong motor and can exert forces of up to 14 
pN (33] • As described by Molineux and coauthors 0, 11^ transcription by RNAP is the 
major mechanism for DNA injection from wild-type T7 into E. Coli and is an intriguing 
additional active mechanism that is of great interest to treat theoretically as well. The 
calculations presented here call for a more systematic experimental analysis of the extent 
to which proteins bind onto phage DNA as it enters the infected cell. 



In this work we have analyzed various effects of DNA translocation of internal capsid 
pressure and "exterior" (cytoplasmic) binding proteins that can be tested by a variety of 
in vitro experiments involving phage ejection kinetics into synthetic vesicles and through 
membranes formed over holes in planar partitions. In these ways one can separately 
control the capsid pressures (by varying salt concentrations or genome length, for example) 
and the nature and concentration of DNA-binding proteins inside the capsid or on the 
other side of the membrane. In addition, it will be important to examine the role of these 
various mechanisms in determining the kinetics of genome delivery in vivo. 
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